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Quantum wires subject to the combined action of spin-orbit and Zeeman coupling in the presence of s-wave 
pairing potentials (superconducting proximity effect in semiconductors or superfluidity in cold atoms) are one of 
the most promising systems for the developing of topological phases hosting Majorana fermions. Other potential 
platforms can be obtained by applying appropriate transformations to the quantum-wire model. One example 
is the recent proposal by Kjaergaard et al. [Phys. Rev. B 85, 020503(R) (2012)], where an effective spin- 
orbit coupling is obtained in electronic systems subject to magnetic textures after applying a local spin-rotation 
mapping. Here, instead, we employ a time-dependent spin rotation that maps the standard magnetostatic model 
into a non-magnetic one where the spin-orbit coupling axis changes as a function of time. This represents 
a proposal for the development of a topological state of matter driven by external forces. From a practical 
viewpoint, though the scheme avoids the disadvantages of conjugating magnetism and superconductivity, the 
need of a high-frequency driving of spin-orbit coupling represents a technological challenge. We describe the 
basic properties of this Floquet system and show that in finite samples it hosts Floquet Majorana fermions at its 
edges despite the fact that the bulk Floquet quasienergies are gapless and that the Hamiltonian at each instant of 
time is a time-reversal symmetric operator. The exact mapping to the static system allows us to show that the 
localized Floquet Majorana fermions are robust to local perturbations, this result is found to be in agreement 
with numerical simulations fully performed in the time-dependent system. 
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I. INTRODUCTION 

Experimental research in the field of hybrid systems involv- 
ing spin-orbit coupling (SOC) and superconductivity is cur- 
rently booming. 1 4 It has been predicted that systems based on 
topological insulators, 5 semiconducting quantum wires' 6 ^ or 
Bose-Einstein condensates (BECs) in ultracold atoms,- 8 ^! can 
realize the topological superconducting phase (TSP) [which 
was theoretically predicted to appear in spinless p-wave paired 
fermions] In particular the one-dimensional (ID) propos- 
als that hosts unpaired Majorana fermions at its edges can be 
thought of a realization of Read and Green phase in or 
of the Kitaev ID modelP-^In these physical platforms the ac- 
cess to the TSP is achieved by breaking time-reversal (TR) 
symmetry with an additional magnetic field, the correspond- 
ing Zeeman energy must be larger than a critical value. 

The virtue of the proposals using semiconducting spin-orbit 
coupled (SOC) samples in proximity to s-wave superconduc- 
tors, is being practical schemes exploding the fact that an ex- 
plicit p-wave superconductivity is not required to realize the 
TSP: as first reported for superfluid ultracold fermion atoms 
systems. 8 While because of the helical nature of the edge 
states in topological insulators the s-wave BCS pairing plays 
the role of a p-wave in spinless systems, 12 in semiconduct- 
ing (ultracold atoms) systems it is the combination of an ex- 
plicit spin-orbit interaction and the s-wave superconductivity 
(superfluidity/BEC) what gives access to the TSP for Zee- 
man energies above a critical value. In finite-size systems the 
success of any of those proposals would lead to an unprece- 
dented breakthrough: the synthesization of the elusive Majo- 



rana fermions (MFs). Beyond the relevance for fundamen- 
tal physics, such a finding — probably already realized — 
would unveil a new set of technological possibilities: since 
space localized MFs are Ising anyons, their topological prop- 
erties can be profited for a pplic ations in quantum information 
and quantum computation! 11 ! 14 ^ 

In parallel, the growing interest in topological phases 
of matter has led to study them in non-equilibrium condi- 
tions. It has been shown, by using Floquet theory for deal- 
ing with time-periodic Hamiltonians, that topological phases 
can be realized including phases that have no analog in 
static systems.^ For example, there has been several stud- 
ies on graphene subject to electromagnetic radiation in the 
microwave-THz regime.^2lln the case of circular polarized 
radiation topological insulating properties arise that give rise 
to the existence of gapless edge states. Another interesting 
proposal is the Floquet topological insulator in semiconduct- 
ing systems. 25 Moreover, ID photonic bound states at the in- 
terface between two distinct Floquet topological phases have 
been predicte dp2l and detectedP^Recently, in Ref.|9] localized 
Floquet Majorana fermions (FMFs) have been predicted to ap- 
pear in cold-atom quantum wires due to periodical change of 
the chemical potential. Similarly, FMF ID modes were re- 
cently predicted to appear at edges of a cold-atom superfluid 
2D system in which the potential of the optical square lattice 
is periodically modulatedP^l 

Our starting point here is the ID model for a quantum wire 
in which spin-orbit coupling, s-wave pairing potential, and 
Zeeman interaction (all of them static and spatially uniform) 
coexist. 6 We intend to gain insight into alternative systems and 
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mechanisms where a topological phase can emerge by apply- 
ing unitary maps to the original model. This is achieved by 
choosing suitable transformations removing either the spin- 
orbit interaction or the Zeeman coupling from the original 
electronic Hamiltonian Hq^, In the transformed system, while 
an interacting term is absent the remaining ones are modified. 
A physical system described by the transformed Hamiltonian 
can potentially share some topological properties of the orig- 
inal system. 28 The removal of the SOC term brings upon a 
successful example of this scheme, 29 allowing the derivation 
of the spin-orbit free topological platform reported in Ref.[30l 
We note that, for systems in which a magnetic field is required, 
it is also important knowing how to replace its effect with al- 
ternative schemesPf^lHere, inspired by concepts of nuclear 
magnetic resonance,^ we report a novel topological platform 
for Floquet Majorana fermions found by removing the Zee- 
man term in the originally static Hamiltonian: a ID quantum 
wire which is subject to a periodical driven SOC. 

The proposed platform does not require any external mag- 
netic field or proximity to magnetic materials. This avoids the 
difficulties of combining superconductivity and magnetism, 
where the Meissner effect could hinder the development of 
topological phases due to the subcritical magnetic screening 
in the quantum wire- see Eq.|5]). This profit is counterpointed 
by the need of a periodical modulation of the SOC-axis that 
should not degrade the pairing potential. In spite of the tech- 
nical difficulties this could imply, our proposal emerges as an 
original case study in the active field of topological phases in 
driven systems! 27 ! 36 ! 

We start by studying the Floquet quasienergy spectrum for 
an infinite quantum wire, finding a family of excitations form- 
ing a Dirac cone that becomes gapless at small wavenumber k 
only for a critical driving frequency Q. c . We predict the value 
of the critical frequency by mapping the Floquet system to 
the static ID topological superconductor system. The system 
here is ungapped in Floquet quasienergies when including so- 
lutions for all k. However, as pointed out by Kitagawa et al. 
in Ref. 18, topological properties can be present in such con- 
dition. 

We also investigate the mean energies of the Floquet solu- 
tions, i.e., the expectation value of the Hamiltonian averaged 
over one driving period. At small k we find that the family of 
solutions in the quasienergy Dirac cone — closing at a finite 
quasienergy value — produces a Dirac cone in mean energy 
that closes at zero energy for the critical driving frequencyP^I 
Remarkably, in this system the topological phase can be dis- 
tinguished from the trivial phase directly from the mean en- 
ergy of the Floquet states at small k. 

We further focus on a finite piece of wire, obtaining Floquet 
Majorana Fermions localized at its edges. This is shown both 
by using the mapping to the static system and by solving nu- 
merically the time dependent problem. The FMFs appear at 
a finite and well-defined quasienergy. Remarkable, the FMFs 
appear even though the Hamiltonian at each instant of time 
does not violates TR symmetry. However, because of the full 
period evolution, TR symmetry is violated and thus the sys- 
tem satisfies such necessary condition for the appearance of 



unpaired Majorana fermions. We have numerically explored 
periodical SOC excitations that do not violate TR symmetry, 
we find that the quasienergy spectrum appears double degen- 
erated and the FMFs are not observed. 

The finite quasienergy at which the FMFs appear is the 
same as that reported in Ref. 9 for one type of FMFs. We 
point out that, because such quasienergy is proportional to the 
driving frequency, the state has zero mean energy. Recently, 
in Ref. |38] Arimondo et al. suggested, with experimental 
support, that the statistics of the non-equilibrium populations 
of the Floquet states could be determined by Bose-Einstein 
(for bosons) or Fermi-Dirac (for fermions) distributions in the 
mean-energy variable. If that is the case the mean energy 
would play the role of the energy in static systems. Similar 
conclusions are theoretically drawn in some limits for a Flo- 
quet system in contact with a thermal bath. 39 These results 
would imply that there might be some regimes in which the 
FMFs can be the highest mean energy occupied Floquet states 
because the FMFs have zero mean energy which is aligned 
with the chemical potential. 

Finally, by working entirely in the time-dependent system 
we numerically verify the robustness of the obtained FMFs 
by showing that they arise in the presence of static local dis- 
order. This is expected by virtue of our exact mapping to the 
static system of Ref. 6 In the presence of such disorder we 
investigate several types of interfaces involving the Floquet 
topological system. 

The paper is organized as follows. In Sec.|II]we outline the 
derivation of the non-magnetic Floquet platform. In Sec. [Ill] 
we explore its properties for both infinite and finite samples, 
identifying the presence of Floquet Majorana fermions. Fur- 



ther discussion and conclusions are presented in Sec. IV 



II. DERIVATION OF TWO ALTERNATIVE 
TOPOLOGICAL PLATFORMS 

We start by introducing the quantum-wire platform in 
which the interactions are constant as a function of time and 
position. In the absence of superconductivity, the electrons in 
the wire follow the Hamiltonian 



O.e 



ll-n: + 11/ ~ "so 



(la) 



2 

pi a 

^kin = =— , H z = E z o- 3 , H so = T p x cr u (lb) 
2m n 

where m* is the effective electron mass, Ez is the Zeeman en- 
ergy, and a is the SOC strength. The <x, (z = 1 , 2, 3) are Pauli 
matrices operating in spin space, while Ez = \gpsB 3 with B3 
a magnetic field applied along the direction of cr 3 . We have 
chosen the most favorable conditions for the development of 
a topological superconducting phase, where the SOC axis is 
orthogonal to be magnetic field. Physically, if the SOC is 
Rashba type the cr 1 axis would be perpendicular to the wire's 
direction x. 

Due to a proximity effect induced by a nearby s-wave bulk 
superconductor there is a nonzero electron-hole pairing char- 
acterized by the energy gap Ao = |A|. Electrostatic gates con- 
trol the chemical potential p.. The system is described by the 
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Bogoliubov-deGennes (BdG) equation®! 
<H = j ^ {x)H^{x)dx, ¥ = (❖^❖J.l'j., -%) , (2a) 

ft \<[,(i)) = ih^ Wt)) , (2b) 
at 

H Q , e -p A 



Ha = 



A* p- H 0M J ' 



(2c) 



where ^^(jr) is the annihilation operator of electrons with 
up (down) spin at position x and time t. In Eq. (j2cj, the quan- 
tities p, A and A* are to be interpreted as multiplied by the 2D 
identity matrix, l2 X 2- The Hamiltonian for the holes, p - ffo,h> 
is obtained by computing the time reversal of ffo,e> 



#0,h - T Hn fi T - CT 2 Ho e (T2. 
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where, 7~ = -icrzK, is the time-reversal operator with K the 
complex conjugation. H^ n and 7/ so are TR-symmetry preserv- 
ing terms in Hn, e , on the other hand, the sign of the Zeeman 
contribution, Hz, is reversed by the TR operator. 

For simplifying treatment on next sections we have writ- 
ten the BdG-Schroedinger equation in Eq.(|2"a|i operating on 
a ket, 10(f)). Indeed, including the case of a BdG time de- 
pendent Hamiltonian, H'(t), the most natural representation 
of those states is a Nambu spinor, <t(x, f) which is written in 
the basis of ^(x, f) and evolves satisfying the Schroedinger 
equation (([' (t)-\hj^ <t>(x,t) = 0. The Nambu field operator 
allows one to write the fermionic annihilation and creation op- 
erator associated with any solution i>(x, f), this will be useful 
below; in particular the annihilation operator written in the 
Schroedinger and Heisenberg picture are 



ft, 



= j $\x, 



0) ■ *¥{x)dx = <5tf(0), 



ft H (f) = 0)M/'(f,Q) , 
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(4b) 



where f U'(t, 0) is the evolution operator from time to f. 
The Hamiltonian Hq of Eq. (2ai has been studied exten- 



sively. It is known that the system is in the topological super- 
conducting phase for Zeeman energies above a critical value, 
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This can be seen by solving Hq<1> = E<t for an infinite wire 
and calculating the Z2 topological invariant that distinguishes 
both gapped phases.^ A topological phase transition involves 
a closing and reopening of the band gap as a function some 
parameter,^ here, the gap closes when \Ez\ - E c z . We get 
insight into the states involved around the critical condition by 
working at Zeeman energies, Ez - E z + Aez- We focus on the 
solutions for small linear momentum p x = fik, the solutions at 
\k\ su kp are gapped due to the s-wave superconducting pairing. 
Close to E = 0, the solutions organize in two branches which, 
at the leading order in the wavenumber k, read 



El(k) = ± y/iak sin <p ) 2 



+ Ae|, 



with 0o = arctan These solutions are gapless only at the 
critical Zeeman field, forming a Dirac cone 



E%k) = ±aksm4>o + 0(k 2 ) = +ak 



0(k 2 ). (7) 



At k + 0, these branches correspond to states that are the 
symmetric and antisymmetric combination of states unper- 
turbed by spin-orbit coupling, with energies E® 2 (k)\ = 



±A<? 7 
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The unperturbed states are <I>, 



(cosf ,0,sinf ,0) and &° 2 = (o,-sin =f,0,cos & ) ^, 
where the Nambu-space spin axis is taken along <x 3 . 

Importantly, for Aez > (Aez < 0) the solution at k = 
with E° > is <£>j (^ 2 '). Passing from one situation to the 
other would imply an abrupt change of such E° > eigen- 
state that can only be achieved after a band inversion: at the 
critical condition the gap closes and both solutions become 
degenerated. It is not possible then to go from a situation 
to the other with a smooth deformation of the Hamiltonian 
without closing and reopening the gap. As mentioned before, 
this indicates that the two gapped phases are topologically 
distinct. 41 Notice that the topological trivial phase is identified 

with the case Ez < -JA^ + p? because the situation includes 
the vacuum phase as such is characterized by a chemical po- 
tential p < and very large. Therefore gapless states only 
appear appear at terminations of the topological phase (with 

Ez > JAq + p 2 ), these terminations may involve interfaces 
with the trivial phase or with the vacuum. 

This kind of midgap states are well-known to appear in 
Dirac -like equations as a consequence of a sign change in 
the mass term. 42 For the BdG equation, due to particle-hole 
symmetry, for each eigenstate with +E a created with the op- 
erator Of, (taken either in the Heisenberg or the Schroedinger 
picture, see Eq.|4|) there is another at -E a which is created 
with the operator O a . 43 Therefore, if a physical system has an 
eigenstate fulfilling the reality condition, i> a = <!)„, i.e., being 
its own antiparticle, it must appear at E a = 0. Furthermore, 
the zero energy state is protected to local perturbations by the 
energy gap when it is found unpaired: localized in space and 
not overlapping spatially with other Majorana solution. Such 
conditions are fulfilled in the TSP phase in the quantum wire 
static proposals. 



In Sec. Ill we shall see that the periodically-driven non- 
magnetic system we treat here shares a similar topological 
transition. Such transition is observed both in the Floquet 
quasienergy spectrum and the associated mean energies. In 
finite samples Floquet localized states at the edges of the sam- 
ple appear protected to local disorder and with a well defined 
quasienergy which is proportional to the driving frequency. 



A. SOC-free magnetostatic platform 



Consider the joint kinetic and SOC terms in Eqs. ( lai and 
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. A spin-dependent shift in momentum, 
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B. Non-magnetic Floquet platform 



F = 

^so — 

U\ =exp(-ife so jccr 1 ), cancels the SOC term. Since U\ rotates 
the spin an angle 2k &Q x around the en -axis, the Zeeman term — 
constant in the original Hamiltonian — changes its axis in the 
(2,3) plane as a function of position. The transformed Hamil- 
tonian, H\# = (j\Hq#U\, becomes, 29 



H 



l,e 



Px 
2m* 



■ E so + E z [cos (2k so x) <x 3 + sin (2k so x) cr 2 ] . (8) 



Recently, in Ref. [30l Kjaergaard et al. demonstrated that 
a system of electrons following Hamiltonian Hi , e in the pres- 
ence of a s-wave superconducting pairing potential, Eq.|2|, 
can develop a topological superconducting phase. Starting 
from arbitrary magnetic textures — generated by, e.g., an engi- 
neered array of nearby micromagnets — and by choosing the 
local magnetic direction as the spin quantization axis, numeri- 
cal simulations show that the effective SOC arising in the local 
reference frame can be sufficient for the emergence of topo- 
logical phases. This results optimal for the sinusoidal mag- 
netic texture of Hi e , Eq. &8\, mapping exactly to the i^o.e of 



Eq. ( la I, as we already showed by applying the inverse map- 
ping. 

Direct comparison with Ho e indicates that the amplitude 
of the rotating magnetic field must be sufficiently strong to 
fulfill the condition of Eq.(|5]l (the chemical potential in the 
equation must be replaced by yU - E so because the energy shift 
E so is absent in this zero-SOC platform). Therefore, the re- 
sults of Ref. 30 supports the approach we adopt here for ob- 
taining new physical topological platforms based on removing 
some interacting terms of interest from a Hamiltonian known 
to have a topological phase by applying unitary transforma- 
tions. Remarkably, the approach opens a complete new fam- 
ily of systems able to host topological properties. This is 
especially interesting in the case of imperfect mappings, as 
for non-sinusoidal magnetic textures in the absence of SOC,^ 
where topological properties can still be present. 

It is interesting to discuss Eq. ([8]) from the point of view 
of spin dynamics of electrons (Ao = 0). Assuming that the 
wavelength of the magnetic texture {2nk^) is fixed the effec- 
tive SOC is more significant to the spin dynamics the weaker 
the magnetic field becomes, namely, in the non-adiabatic spin 
evolution regime. Obviously, if the magnetic field amplitude 
Ez is fixed the same effect is achieved by decreasing the length 
2nk~Q . In these ways, spin eigenstates of Eq. ^ are not fully 
aligned with the local magnetic field, otherwise no spin mix- 
ing survives from the position rotating magnetic field. The 
adiabatic regime is defined in the limit oj s » 2n/t c , where 
u) s = E z /h is the Larmor frequency of spin precession and 
f c = L/vp is the time it takes a spin carrier with Fermi velocity 
vp to cover a characteristic length scale L over which the di- 
rection of the magnetic field suffers a significant change (see 
Ref.l45land Ref.|46]for a detailed discussion). Non-adiabatic 
dynamics requires w s ~ 2n/t c . 



Here we get rid of the Zeeman term of Eq. ( lai by apply- 
ing a global time-dependent rotation to the spins along the 
magnetic -field axis with an appropriate frequency. For sim- 
plicity we first consider the electron block of the Hamiltonian 
(not the BdG equation) and thus the mentioned transformation 
is U r (t) = exp {-icr^tEzIn). The time evolution of the rotated 
states 

|0 rot (f)) = f/ r f (f)|0(O) , (9) 
given by the equation \hj t \(f> mt (t))-H e (t) |0 rot (f)), with 
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Where we have used that the equation of motion for the origi- 
nal states, ihj t \<p(t)) = #o,e(0 \<f>(t))- Notice that states \(f> Iot (t)) 
are, indeed, the interaction-picture states resulting from con- 
sidering H^ n +H so as a perturbation to flz- This means that 
a physical system following the time-dependent H e (t) in the 
Schroedinger picture can be mapped into the static ffo,e by 
the unitary transformation U r (t). Since topological phases can 
arise from Ho e when a pairing potential is present, we ex- 
pect similar features from the periodically driven H e (t) in the 
framework of Floquet theory. 

It is clear that the excitations in systems described by ei- 
ther H e {t) or Ho,e are not identical when each of them is as- 
sumed to be in the Schroedinger picture. (Even in absence 
of spin-orbit coupling, where the H e (t) results to be time in- 
dependent, actually both Hamiltonians have different energy 
spectra). However, assuming electron-hole s-wave pairing, 
both systems share the topological trivial phase. Similarly, 
for the non-magnetic system with a periodically rotating spin- 
orbit axis, our mapping allows to predict the existence of a 
topological phase. As we shall show, the critical magnetic 
field in the original Ho, e implies the existence of a critical 
driven frequency in the Floquet system above which a finite 
piece of wire hosts two Floquet Majorana states localized at 
its edges. Furthermore, these edge Majorana solutions appear 
even in the presence of disorder despite the fact that the Ma- 
jorana Floquet states are not gapped in the quasienergy spec- 
trum. 

The fact that states responding to the Hamiltonian H e (t) can 
be transformed by a time-dependent unitary mapping to the 
states responding to H(, e is not enough to derive topological 
properties in HJf), Here, the reason why topological prop- 
erties are also expected in the time-dependent system lies in 
the simple nature of U r (t). In general, any BdG Hamiltonian 
H@(t) — including those without topological properties — can 
be mapped into the topological one Hq of Eq.((2]i (the one 
based on ffo.e)- This is possible provided the unitary trans- 
formation is sufficiently complicated to the point of introduc- 
ing the requirements of the topological phase. More explic- 
itly, starting from the solutions to the arbitrary Hamiltonian, 
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(H@{t) - ifij t ) |0 e (O), the topological Hamiltonian Ho gov- 
erns the time dependence of the transformed states if and only 
if 1^(0) = <U top (t,0)<Ul(t,0)\<f> & (t)), where 0) and 

f Ui QV (t,Q) are the time-evolution operators associated to H& 
and Hq, respectively. In our case, Hq = H(t) [which is the 
BdG Hamiltonian based on H e (t) in presence of a pairing po- 
tential, see Eq.(fT2]i], and 



U K u) = K uv uA)YUl(t,0) = [ 
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(11) 



is a rotating-frame spin transformation that does not introduce 
topological superconducting features. Therefore, topological 
properties in both systems, if any, must be linked. In the fol- 
lowing section we show that the solutions of the non-magnetic 
BdG Hamiltonian based on H e (f) in the Schroedinger picture 
present topological features related to those of Hq. This is 
because U^(t) establishes a simple mapping between the en- 
ergy spectrum of the static Ho and the Floquet quasienergies 
of H(t). Furthermore, we shall see that each eigenstate of Ho 
is associated with one Floquet quasi-energy state sharing its 
properties. 



III. PROPERTIES OF NON-MAGNETIC PERIODICALLY 
DRIVEN SYSTEM 

We consider a system of electrons described by the time- 
dependent Hamiltonian of Eq. ( fT0) i in the Schroedinger pic- 
ture. The rotating SOC coexists with a constant s-wave pair- 
ing potential, where the corresponding BdG equation reads 



(x)H{tyV(x)dx, 



, /m _,#e(0-j" A \ 



(12a) 
(12b) 



For convenience in the following, we choose the supercon- 
ducting phase equal to zero, A = Ao, we introduce the driving 
frequency, Q., a time-dependent SOC vector, A(f), and a static 
SOC vector, A (0) , rewriting the electron and hole Hamiltoni- 
ans as 



H e (t) = H h (t) 



Px 

2m* 



+ ^(A(0 + A"»). < r, (13) 



where er = (cr\ , <7"2, 1x3). Equation ( fT3j ) follows from Eq.Q af- 
ter noticing that, for any given time t', H e (t') is a time-reversal 
symmetric operator. For the ideal rotating (IR) case — which 
is the best candidate for hosting Floquet Majorana fermions — 
the SOC vectors are, 

A(f) = Am(0 = (a cos (Q.t) , -a sin (Qf) , 0) , (14a) 
A (0) = = (0,0,0). (14b) 



Notwithstanding, we shall also discuss some scenarios away 
from this ideal situation. 



Periodically driven Schroedinger equation and Floquet systems 



As discussed in Sec II B a ID system which is subject to 
a rotating SOC is related to the magnetostatic, topological 
quantum wire. Time-dependent periodic systems are better 
described in the context of the Floquet theory. Here we briefly 
describe such theory. The starting point is the need to solve 
the Schroedinger equation, 



(15) 



H(t)-in-\m)) = o, 



for a Hamiltonian periodic in t, H(t) — H(t + T). Energy is 
not conserved, however, by virtue of the Floquet theorem the 
solutions can be written as 



\Ut)) = z-* e * t \<t>l{t)), 



(16) 



where T = || is the driving period and |<^„(0) is a peri- 
odic state and the subindex a encodes the quantum numbers 
of the different solutions. We define the Floquet operator 
H F = (H(t) - itij,)- As \(p a {t)} is a solution of the Schroedinger 
equation then it must hold H F \<p a {t)) = 0, thus one finds 



/-/in - m- 1 \cf> T a (t)) = s a \<t> T a (t)). (17) 



The quasienergies, e fl , are the eigenvalues of the Floquet 
operator Hp. The corresponding eigenvectors are the Flo- 
quet quasienergy states (QESs), ^^(t)^. Besides, the evo- 
lution operator associated to H(t) has the form *ZY (t , fo) = 
T f exp^-|J^ H(t')dt'^, where T, stands for time-ordering 
symbol. From Eq.(fT7]i it then follows 



<U {to + T, to) \<t> T a {to)) = e-* e " T \<p T a (t )) . 



(18) 



Hence, the quasienergies can also be extracted from the phase 
factors er'r> e ° T : eigenvalues of the evolution operator over a 
driving period. 

Floquet theorem is equivalent to Bloch theorem for pe- 
riodic driving instead of position periodic potentials, where 
quasienergies in Floquet QESs play the role of momenta in 
Bloch states. The equivalence of solutions belonging to dif- 
ferent Brillouin zones in momentum is one of the main signa- 
tures of Bloch systems. Similar properties are shared by Flo- 
quet systems in the quasienergy axis. This can be seen from 



Eq. 16 after noticing that states \(f> a (t)) remains unchanged by 
the substitution 

s a ^s a +nm, \<[> T a (tj) -> e i,,n ' \<f> T a (tj) . (19) 

Namely, the new Floquet state |0a'"~ shlft (f)) has shifted 
Fourier components so that for each physical state there 
are infinite Floquet QESs that are nonorthogonal, with 
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Ul(t)\f' n ~ sim (tj) = e'""'. However, orthogonality is recov- 
ered by defining the inner product of Floquet states as the stan- 
dard inner product averaged over a driving period, 



((timid))) = I JT {fiit'Mit')) 



dt' = 5 a j 



(20) 



To retain only one QES for each physical state it is sufficient 
to choose any range of quasienergies of length HQ. — such as 
[e, e+HQ.) — and discard all the QES outside this region. 
The time-dependent Eq. ( 17 1 for the Floquet states still 



needs to be solved. One possibility is to compute the evo- 



lution operator and solve Eq. ( 18 1. We do this for the case of 



a finite piece of wire. Alternatively, we can switch to Fourier 
representation by rewriting the Hamiltonian and the Floquet 
states as 

oo oo 

H{t) = ^ e -'"'"#"> , \f a (t)) = J] e-' fi "'|^). (21) 

n=— oo n=— oo 

By direct substitution in Eq. ( [T7| ) we find 
£,„ (flW - nttMoJ \<fft~ m) ) = e a \<f,^), which in ma- 
trix representation reads 



H m H (0) H (-D 



ucoy 



(22) 



This configures a time-independent infinite dimensional 
eigenvalue problem where s a and are unknown. When- 
ever possible, we solve this analytically provided some ap- 
proximations are introduced, by comparing HQ. with other 
energy scales appearing in the Fourier components of 



Eq. (22 1. We shall do so only in a small (but important) re- 
gion in parameter space. For dealing with general situations 
we resort to a more exact numerical approach. 



Mean energy of Floquet quasienergy states 

The mean energy of a Floquet QES is defined as the ex- 
pectation value of the Hamiltonian averaged over one driving 
period: 



- i r T 

E a = - J dt(<f> T a (t)\H(t)\<p T a (t)) 



T 

— s„ — Q 



dtU' a (t) 
'<9Q' 



d d 
H(t)-ih—+ih — 
dt dt 



(23) 



where we use an extension of the Hellmann-Feynman theorem 
(for the Floquet operator Hp)W^ Notice that computing the set 
of quasienergies s a by fixing Q. — the usual procedure in most 
numerical treatments of Floquet systems — does not permit 
the calculation of the associated E a . To this aim one needs to 
either know ds a /dQ. or use the periodic time-dependent eigen- 
state \<p\ (f)) to compute the integral given in Eq. (23 1. 



Starting from a Floquet state |<^J(f)) with quasienergy e a , 
we find that the mean energy is identical for all associated 



shifted states of Eq. ( 19 1 with quasienergy s a + nftQ, since 



E n a shl " = (s a + iitel) -n-^-(s a + ntiQ) = E a . 
oil 



(24) 



In summary, the mean energy is a useful quantity for classify- 
ing the Floquet QESs and the physical state \<p a {t)) of Eq. ( 16 1. 



Furthermore, very recently, in Ref. [38] Arimondo et al. elabo- 
rated on the hypothesis — while comparing with experimental 
results with apparent success — that the population of Floquet 
states might follow Bose-Einstein (for bosons) or Fermi -Dirac 
(for fermions) distributions in the mean-energy variable. The 
latter is an additional motivation towards the study of the mean 
energies in Floquet systems. An important point we shall ad- 
dress here is how the mean energies of Floquet states behave 
at both sides of a topological transition. 



A. Bulk properties 

We now investigate the Floquet solutions of the time- 
dependent system in its infinitely long version. For arbitrary 
periodic driving, one must proceed numerically because com- 
plexity impedes analytical treatments. For the ideal rotating 
case our numerical implementation is in full agreement with 
the analytical results that, as we show in detail below, follows 
from the mapping to the static topological problem discussed 
in SecHTBl 

As the system is translational invariant the linear momen- 
tum, p x = Hk, is a good quantum number; in the follow- 
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(a) E sn = (b) £,„ = 0.25A (c) E m = 1 .75A„ (d) Ej— 



ing we use the wavenumber k to label the solutions. For a 
given k, Hkit) is a 4 x 4 matrix including electron-hole and 
spin degrees of freedom: the Nambu space. The trial state 
\4>k{t)) = \k) \<Pk(t)) m the Schroedinger equation ([15} leads to 
a differential equation for the Nambu ket, \<fik(t))'. 



(25) 



B&)~in-Mty> = H$ i^(t)> = o. 



For each value of k there exist four solutions, we write them as 
|^/t,a(0) with a = 1, ... ,4. The knowledge of \(fik, a (t)) allows 
us to write the solutions as in Eq.(fT6| by virtue of the Floquet 
theorem for the periodically driven Hk(t). 

In order to compute the Floquet solutions we extend the 
numerical method for the calculation of dispersion relations 
in spatial lattices given by Ando in Ref.|49] This is performed 
by discretizing time along a driving period. For each tem- 
poral "site" i located at f,- = idt the local Hamiltonian //*(*,•) 
acts on the Nambu spinor ^(f,). The discretized version of 
the term — ihj in Eq.(f25]l provides a hopping term between 
first-neighbors sites. The discretized //* can be thought as an 
effective spatial (in the "sites" f,) periodic lattice Hamiltonian 
in Nambu space. The momentum associated to the motion 
on the lattice f, is the Floquet quasienergy, s. Therefore, the 
Bloch solutions in this effective lattice are the Floquet ones 
in the time-dependent system, so that by imposing Eq.(f25)l 
one obtains e as the variable that would be the momentum in 
a spatially periodic lattice. The technical details of the ex- 
tended Ando method for Floquet systems will be presented 
elsewhere. The results are equivalent to those found by com- 
puting the evolution operator over a period of the driving fre- 
quency and then solving the eigenvalue problem of Eq.(18). 



1. The IR driving: Floquet bands and Majorana states. 

Figure [T] shows numerical results for the time-dependent 
SOC of Eq.([14), where the SOC axis rotates harmonically 
around <x 3 axis. As a starting point, in Figs. [TjT.a), [TJC.a) 
and[TjB.a), we show the dispersion in the absence of driving, 
corresponding to the eigenenergies E(k) of the static part of 



the Hamiltonian tr® of Eq. (21 1. The only interacting term 



in H m is the pairing potential, so that electron-hole states are 
spin-degenerated with eigenenergies 



E ±M = ± ^Al + (h 2 k 2 /(2m*)-fi) 2 



(26) 



These appear as solid lines. Additionally, the same disper- 
sions shifted by an integer multiple of TiQ. are displayed as 
dashed lines. Those bands are relevant for the time-dependent 
system since, as shown in Eq.(p2]i, the W- m > operators (m^O) 
mix solutions of i7 (0) with energies differing in mtiCl. 

The Fourier decomposition of the interacting term p A -Aj R (f)- 
<r/h appearing in both H e (t) and H^t), with Air(?) given in 
Eq.(fT4|i, determines that the only off-diagonal terms contribut- 
ing to Eq. (|22]> are 



rV(=Fl) _ rV(=Fl) _ & „ _ 

H e -H h - j i PxCr ± , 




-3 3 -3 3 -3 3 



(C.a) 



(C.b) |(Q Q.-Q.J ( C . C ) 



(C.d) 




ilk 



ilk. 



iie 



(27) 



FIG. 1. Ideal rotating (IR) case of Eq.(|l4}: Floquet quasienergies 
and mean energies. Reference energy and momentum are Ao and fc A = 
V2m*Ao/S, respectively. We set yu = 2Ao so that 7ifl c = 2 V5Ao. Top 
(T), center (C) and bottom (B) panels show results for CI = 0.9O c , 
fl = £l c and £1 = l.lfl c , respectively. Dotted horizontal lines at 
nKl are included in all panels. Panels (a) show the spin-degenerated 
naked dispersion E ± k<T of Eq. ( |26} in the absence of the IR SOC-field 
(solid lines). Dashed lines correspond to E+ ^+mnCl branches which 
are important for the mixing due to the periodic driving, see Eq.(|22j. 
(b) and (c) Floquet quasienergies s for finite IR SOC-field of inten- 
sity E so = a 2 m*/(2h 2 ) = 0.25A and 1.75A , respectively. Numerical 
results show an excellent agreement with analytical calculations ob- 
tained by mapping the time-dependent problem to the static system 
having a constant SOC and a finite Zeeman field: the solid black lines 
depict the four bands of the corresponding static system (shifted in 
energy by +n£l/2, see text). Panels (d) show the mean energies E 
of the Floquet states for the simulated IR SOC-field amplitudes. In 
(C.b) and (C.c), two branches of Floquet solutions close in a Dirac- 
like quasienergy cone for small k. The reopening of the Dirac cone 
for CI > Cl c , (B.b) and (B.c), indicates a topological phase distinct 
from that for CI < Ci c . The mean energy shows a distinct pattern on 
each phase (see text and Fig|2j. 



where cr ± — <j\ + io"2 are the <r 3 -spin raising and lowering 
operators. This means that the driving introduces mixing be- 
tween H (()) +mfiCl blocks with (m-m') = 1. Due to translation 
invariance, for each k one needs to solve a time-independent 
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matrix Floquet equation as Eq.(|22|>, where each block is a 4x4 
matrix. The amplitude of the IR mixing terms is +ak. This 
means that at k — the system is static and the Floquet treat- 
ment is not required. However, in the neighborhood of k = 
the system does need Floquet treatment. The energy of the 

solutions for k — > are ± ^Ja^ + jj 2 . For a ftQ. = 2 ^Ja^ + fi 2 , 
the energies of Floquet blocks with (m-m') = 1 cross at k = 0. 
In particular, the (m, -) band crosses the (m— 1, +) one. This 
situation is illustrated in Fig[TJC.a) as the upper band in solid 
black line (m? = 0) crosses the lower band in dashed black line 
(m = 1). From this condition we define the critical frequency, 
Q. c , as, 



Q, 



(28) 



This £2 C , determined here within Floquet theory, coincides 
with value predicted by the mapping to the static case in- 
troduced in Se c|IIB| [starting from Eq.([T0|> by identifying 
Q=2£ z /ft,Eq.(|^pollowsfromEq.((5J]. 

In Figs. [TJCb) and [TJC.c) we show the resulting Floquet 
quasienergies at the critical frequency for two different ampli- 
tudes of the rotating SOC. For small k we obtain Dirac-cone 
like quasienergy excitations closing at values 



£ 0> „(Q) = (n + 1/2)M2, 



(29) 



with = Q c . They coexist with other bands also closing at 
£o,«(ficX but with vanishing ds/dk. The bands do not longer 
close at k » for either Q < Q c [shown in Fig[TjT.b) and 
FigJTjT.c)] or Q > Q c [shown in Figs. [TjB.b) and [TjB.c)]. 
Such opening and closing of quasienergy bands suggest a 
change in topology. However, here the full quasienergy spec- 
trum is gapless as the solutions for larger k take all possible 
quasienergy values. 

At this point it is clarifying to make use of the mapping 



introduced in Sec II B where the Floquet solutions can be ob- 
tained from the static solutions of the SOC plus Zeeman ID 
superconducting system. We start from the eigenstates of Hq 
with energy E a [see Eqs. ( lb i and (2a i] 



WW) 



(30) 



with S> a (x) a Nambu spinor satisfying (Hq - E a )iP(x) = 0. 
The Hq eigenenergies for the infinite case are distributed in 
four bandsP 



E a = E bub2 (k) = bi 



A 2 + E 1 , - A + <rV 



2m 



(31) 



1/2 



with b\ — ± and = ±- The spin rotation of Eq.([TT| applied 
on the static solutions give 



(32) 



which are the solutions of the non-magnetic IR SOC su- 
perconducting ID system by substituting Ez — > fttl/2 and 
a — > a. 

The physical state \(p a (f))< Eq.([16]l, can be associated to a 
family of Floquet QESs, the members of such family have 
quasienergies differing in multiples of ftQ. Disregarding 
such an ambiguity, Eq.(|32]> is important to determine the 
quasienergy values associated to the solution with energy E a 
of the static quantum-wire. For convenience, we first write 
each eigenstate of the static system of Ref. 6 (assumed known) 
as a Nambu spinor, <£>"(x) = (A a (x), B a (x), C a (x), D a (x)) T . 
We cast the resulting physical state in the form \(f> a (t)) = 
e ~T' e "' |</>J(f)), this allow us to identify of the quasienergy for 
the time-dependent system as 



s a = E a + (n + l/2)»fl = E a + e ,„(Q), 
and the associated Floquet QES as the Nambu spinor 



(33) 



AnClt 



( A a (x)e nt \ 

B„(x) 
C a (x)e in < 
D a (x) 



(34) 



Importantly, Eqs.([34[i and ( |33[ l are valid also for cases in 
which both linked systems are not translational invariant and 
thus the solutions are not eigenstates of the momentum p x ; 



for example in Sec III B we deal with finite and disordered 
systems. 



We set n — in Eq.([33[> to generate Floquet quasienergies 
of the IR case using the energies of the static system given in 
Eq.((3T|. We plot in black solid lines the resulting s a in pan- 
els (b) and (c) of FigfTJ i.e., for different driving frequencies 
and Rashba strengths. Floquet quasienergies obtained entirely 
within the Floquet picture with our numerical method show an 
excellent agreement with the analytical results. In Figs.[TJC.b) 
and[T|C.c) for Q = £2 C we see that the Floquet-Dirac branches 
closing at sq „(Q. c ) map to the Dirac-cone branches of Hq that 
close at E a - (those with nonvanishing ds/dk at k = 0). In 
the plot, we have chosen /j > Aq to show that the mapping 
holds in general. It is known that the /i = case is the most 
favorable for expressing the Kitaev model in the static Hq be- 
cause the effective p-wave gap is largerFS Here, y. > Aq as- 
sures a nonzero electron density allowing the development of 
superconductivity by proximity effect when switching off de 
driving (a-0). 

In the static system, as discussed above, the closing and 
reopening of the band gap indicates a change in topologi- 
cal properties. There, when MFs are present, they appear 
as solutions of the BdG equation with E a = Emf = 0. This 
follows from the particle-hole symmetry of the BdG Hamil- 
tonian that enforces zero energy for any eigenstate which 
is its own antiparticle. In the translational invariant system 
E a = solutions only appear at the critical Zeeman energy 
for k = 0. Either in extended or finite systems, the general 



form of a Majorana fermion in Nambu space is (?'^ [F (x, t) — 

(^mf(*X Bmf(jc), ^mf( x )' "^mfW) ' an( ^ tmis assoc i ate d 
fermionic operator is (because E a - both the Schroedinger 
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and the Heisenberg picture lead to the same operator, see 



6° = 6° 



+B MF (x)V l (x) + B* MF (xWl(x)] . (35) 



A natural question arising here is whether the existence of 
Majorana solutions in the static system dictates the existence 
of Majorana solutions in the IR system of Eq.( |12a| >. To see 
this, we construct the solution of the driven IR system which 
is associated to the static MF solution of Ho; in the Nambu 
spinor representation this is just <P F mf(x, t) - U^(t)0 ( l AP (x, t), 
we get 



<Z>fmf(*, 



A MP (x)e 




B MP (x)e- 




B* MF (x)e 






-Lot 
2 I 



(36) 



we see that this state satisfies the particle-antiparticle condi- 
tion while it evolves periodically in time [this follows from 

Eq.Q noting that & H (i) = <G'\t, 0)O H (0)'ZY'(f, 0)) = & H (ff]: 
it is a Floquet Majorana fermion (FMF). Notice that two driv- 
ing periods are needed for the time-dependent state to revisit 
the same instantaneous configuration. This is so because the 
t/R(f) of Eq. ( 11 1 is made of spin- 1/2 operators with rotation 
angle 2nt/T (see Sec III A 2 1. For t = 2nT (rotation angle Ami), 
it reduces to the identity. This is consistent with the fact that 
a representation within Floquet theory of such a MF state — 
E a = in Eq.([33]> — must have quasienergy 



sfmf = £o,n(m = (" + 1/2)7*0, 



(37) 



meaning that FMF states acquire a phase of n+2nn (i.e., a fac- 
tor -1) after one driving period. At £2 = Q C , the latter coincide 
with the quasienergies at which the Floquet-band gaps (those 
with k^O) close. 

Similar FMFs with finite quasienergy were reported by 
Jiang et al. in Ref. [9] where they studied a BEC cold-atom 
quantum wire subject to static SOC and magnetic field with 
a periodically modulated chemical potential. Their results 
demonstrates that finite-quasienergy FMFs are not restricted 
to our particular model. An important point that Jiang et 
al. brought upon is that the BdG-Floquet operator preserves 
electron-hole symmetry. In his context this means that for 
each Floquet state with quasienergy +e a created with the op- 
erator & a (see Eq.Q) there is another at — s a which is cre- 
ated with the operator <J> fl . Because quasienergies differing in 
nWl lead to equivalent physical solutions this means that the 
Majorana condition (Oj, = i> fl ) can also be satisfied by states 
at exactly the nonzero quasienergies given in Eq.((37]>. We 



shall return to the discussion of unpaired FMFs in Sec IIIB 
as bounded states appearing at the edges of topological sys- 
tems are intrinsic to the finite size samples treated there. 



2. The IR driving: topological properties 

We have successfully established a link between the 
eigenenergies E a (eigenstates !*„(?))) of the static quantum 
wire of Ref. 6 , Ho, and the quasienergies s„ (Floquet states 
\</>a (0/) °f our proposed non-magnetic system, H(f). We have 
shown, using the exact mapping, that for each Majorana so- 
lution appearing in the static system there is an analog Flo- 
quet Majorana fermion appearing in the driven system. Those 
FMFs are protected by the same energy gap that the associated 
MFs in the static system. 

In order to make even more explicit that the topological 
properties of the Floquet system are equivalent to those in the 
static system we make use of our mapping to write the evolu- 
tion operator for the rotating states. From Eq.((32]> at t — the 
initial state in both pictures is the same, 



\<pm) = |*2(0)) ■ 



(38) 



Because of the time-independence of Hq we can write the time 
evolved state in the static system as, 



|*°(0) = exp (-iA r) |02(O)). 



(39) 



In order to simulate the time-dependent system with the IR 
excitation, the amplitudes of the interactions in Eq.(fTb"]i are 
chosen as a — > a and Ez — > Ml/2. Introducing Eq.(|3~8]) 
and Eq.([39]> in Eq.([32]> we write the evolution in the time- 
dependent system as, 



l* B (0>=W,Q)|* B (0)>, 



(40a) 
(40b) 



After replacing Ez — > n£l/2 in Eq.(fTT]i we get (/^(T) = -I4, 
this is a 2n angle spin rotation, with I4 the identity in the 
Nambu space (in the same basis as Hq given in Eq.pcj)). 
Therefore, when particularizing f U(t,Qi) for the evolution over 
a full period we get 



(41a) 
(41b) 



'Z/(7\0) = exp|--H eff r 
H eff = H + tiSoJQ), 



This means that the evolution over a period in the time- 
dependent system, H (T, 0) = T, exp (— ^ J H(t')dt'), is equiv- 
alent to the evolution over the same period T with the time- 
independent effective Hamiltonian, H e g . Not surprisingly, the 
effective Hamiltonian is just Ho trivially shifted in energy by 
the amount so M (Q.) = (n+ l/2)hQ. defined in Eq.(|29]i. In the 
light of Eq .(ff8] l, such shift justifies the form of the quasiener- 
gies in Eq.(|33|>. 

As pointed out by Kitagawa et al. in Ref. 18 the topological 
properties contained in the effective Hamiltonian of a Floquet 
system are to be studied with the tools for time-independent 
topological phase transitions. The topological phases appear- 
ing in H e ff imply associated topological phases in the Floquet 
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system. We can now formally conclude that the same topo- 
logical invariant, a Z 2 number,^ which is different in the two 
gapped phases of the static and magnetic model of Ref. is 
also identifying the two Floquet topological phases in the non- 
magnetic driven system. In particular this corroborates that 
the regime with > fi c (supercritical) is indeed a Floquet 
topological superconducting phase whereas for CI < Q c (sub- 
critical) is to be associated with the topological trivial super- 
conducting phase. 



3. The IR driving: mean energies 

We now discuss the mean energies of the Floquet states de- 
picted in Figs. [TjT.d), [TJCd) and[TjB.d). These follow from 
Eq.|23) after calculating the derivatives de a /dQ.. The numer- 
ical results agree with the analytical mean energies derived 
from Eq.pTI by computing 



Efl En E\ 



dE a 
: dE z 



(42) 



and substituting Ez — > Ml/2. At critical frequency, Q. - 
Q. c , the two Floquet states with quasienergy forming a Dirac 
cone — see black solid lines in Figs.[TJC.b) and[TJC.c) — lead 
to linear mean-energy dispersions around k — as shown in 
Fig. [TJC.d). We verified that this happens as a general rule, 
where the mean-energy bands closing the gap at E = derive 
from the Floquet-quasienergy bands closing at eo,«(^c)- 

As a result of a discontinuity in the mean-energy bands as 
a function of the driving frequency, these get very sensitive at 
fi c for small k. Such a discontinuity is intimately related to 
the change of topological properties. Right at k-Q, the ampli- 
tude of periodic driving vanishes and the mean energy simply 

reduces to the energy, E ± (k-Q) = ± ^Ja^ + /j 2 = hQ. c /2. Below 
the critical frequency, the top mean-energy band is connected 
to the E + (k = 0) point while the bottom band acts likewise 
with the E-(k = 0) point, instead. Remarkably, above the criti- 
cal frequency this relationship is inverted. The inversion point 
is defined at Q c , where linear dispersion in mean energy is 
developed excluding the states with k-Q: points E ± (k = 0) 
appear isolated from the bands (see following discussion and 
Fig. [2] for details). In general, for higher k the mean-energy 
spectrum is complicated, this is particular for the case shown 
fi > Aq. For n « the mean-energy bands very much resem- 
ble those of the energy of the static system for Q. = Q c , away 
from the critical condition the low k mean energies are always 
different than the energies of the static systemP3 

For a better understanding of criticality in the vicinity of 
Q c , it is useful to solve the Eq.(|2"2"|) for the Floquet oper- 
ator in the limits of // « Ao and small k. We first need 
the solutions of the static component of H{t), i.e., i/ (0) , 
in the Nambu basis with spin quantized along the 3-axis, 
{xfr^, ij/^, ifri, —ifrffi. For simplicity, we define Nambu unit vec- 
tors as tf/[ = {811,821, Sjj,84j) T . As 7/ (0) only includes the 
kinetic term and the pairing potential, we define the angle 



sarctan(Ao/(ft 2 £ 2 /(2m*) - fi)). The solutions are 



l+,*T> 



VZ 

^'ikx 



fcos + sin —tlit 



<Pk-_ 
2 



2 



) 

e"" / <fk— <Pk-~\ 
+, k J,) = — I cos —\]/ 2 + sin — ij/ 4 1 



,*T> 



VZ 

e ikx 

n 

Jikx 



^sin - cos 



<Pk-. 
2 



<A 3 ) 

, .. e"" / . <p k ~ tp kr \ 
VZ ^ m 2~ ~ 2~ ' 



(43a) 
(43b) 
(43c) 
(43d) 



with the energies given in Eq.(|26]>. These states are, indeed, 
eigenstates of each block H^+mMl of the Floquet operator 
given in Eq.p2]). We define \m, +, kcr) = e™ n ' |±, kcr) for re- 
ferring to the eigenstates of block m by simply introducing a 
Fourier phase factor. 

We limit our analysis only to blocks m = and m = 1 in 
Eq.(p2]i, assuming the mixing of few bands only. This is jus- 
tified since: (i) the only mixing components in H(t) due to 
the periodic driving are i/ (±1) , (ii) in the small k limit the 
time dependent components are small, discouraging mixing of 
states far apart in energy, (iii) close to the critical frequency, 
Q. = D. c +6e/fi, other states with small k are at least Ml c away 
in quasienergy, and (iv) other choices with m—m' = 1 lead to 
shifted quasienergies (in multiples of Ml) and, therefore, to 
equivalent solutions. Furthermore, we restrict the Floquet op- 
erator to the subspace of |0, +, kcr) and 1 1, -, kcr), i.e., the four 
states crossing when Q = Q C . 

Combining Eqs. |27| and ( p~2] > up to linear order in k, the 
Floquet operator of Eq.d22| reduces to 



[ no.j2 000 

M2 c /2 a£sin<£>o 

ak sin ipo M2 c /2 + 5e 

£Q c /2 + 5e ) 



(44) 



in the ordered basis {|0, +, k T> , |0, +, k I) , |1, -, k f> , |1, - k J,}}. 
We readily notice that the 1 st and 4th states 



|0[(f)) = \+,k t> with e, = ^Q c /2, 
\<Pl (0) = 



(45a) 



k i) e 1 "' with e 2 = MIJ2 + 8e, (45b) 



do not mix under the action of the rotating SOC axis. Their 
mean energy, Eq.(|23|>, is identical to their energy: 



(46) 



Ei = hQJ2 , E 2 = -Ml c /2, 



where we used that dQ. c /d£l = and d(6e)/dD. = h. The 
remaining two solutions, for k + 0, are a mixture of states with 
different m. By defining the angle % = arctan {2ak sin (po/Se), 
these solutions 



\$(t)) = 



cos^|+,A:|>-sin|K A;T>e i£lf , (47a) 
sm^|+,U> + cos^f|-,JtT>e i0f , ( 47b > 



11 



with quasienergies, 



1 / 6 ^ 

s 3 (k) = h£lJ2 + -Se - + (aksimpo) 2 , (48a) 

^(y) 2 + (^sin^ ) 2 . (48b) 



e 4 (k) = m.J2 + -6e + 



At the critical frequency, 5e - 0, the two branches form a Dirac 
cone in quasienergy. The vertex at k — lies at quasienergy 
SQ c /2 (while the /-shifted solutions cross at hQ. c /2 + lTiQ. c = 
e/,o(fi c ))- The superconducting gap is essential since the mix- 
ing term aksinipo, with sini/?o = ~Ao/ + fi 2 , vanishes for 
Ao =0 (leading to fiat dispersions). 

For arbitrary values of k (always in the small k approxima- 
tion) the mean energies for these two solutions are 



E 3 (k) = ~E A (k) = 



6eh£l c /2 - 2 (ak sin <po) 
^5e 2 + 4 (ak sin <po) 2 



(49) 



At the critical frequency, 6e - 0, the two branches form a Dirac 
cone in mean energy, +\ak sin <p$\. This excludes the case k 
0. with mean energies fixed at +hQ. c /2. We find 



E 3 (0) = -E 4 (0) 




for Q < Q c 
for Q > £1, 



(50) 



In Fig|2]we show the different behavior of the quasienergies 
and mean energies involved around the transition, i.e., the so- 
lutions 3 and 4 in the latter derivation. For simplicity only the 
bottom quasienergy [mean-energy] band is plotted in the plane 
(Q, k) in Figj2ja) [FigQg)]. The low-A; quasienergy band 
passes through the transition without visible change (see up- 
per panels in Figj2jb) to Figj2jf)); on the other hand the mean- 
energy band (lower panels in Figj2jb) to Fig[2jf)) clearly al- 
lows to distinguish the two sides of the transition. Both the 
abrupt change in mean energy and the Dirac -cone dispersion 
at the critical parameter are due to the band inversion and 
linear-in-fc mixing underlying the phase transition. The re- 
sult here is consistent with what is commonly found in Flo- 
quet systems, i.e., avoided crossings in quasienergy producing 
abrupt changes in the associated mean energiesPH 

In our study of the Floquet problem, the need to compute 
the mean energies E a arises naturally. The emergence of the 
mean energy as an indicator of the band inversion is remark- 
able. It is therefore pertinent to ask whether E a is a mean- 
ingful quantity to be computed in static systems those dis- 
cussed in Ref. |6]or even in more general contexts. We have 
shown that E a can be derives from Eq.(|42]> as a function of the 
eigenvalues of the static system. By choosing an eigenstate of 
Ho, |</>fl(?))i and applying the Hellmann-Feynman theorem we 

find that E a is just the expectation value in state (^"(f)) of the 
Hamiltonian operator without the part proportional to E z , i.e. 
without the Zeeman interaction. In other words E a is the ex- 
pectation value in |0a(f)) (eigenstate of the full Hamiltonian) 



£'=e-(n+l/2)h£l 



(£\*,£i) =(o,o, a) 




( b)0.8Q c ( c)0.995fl c (d)fl c ( e)1.005fl c (f)1.2fl c 
^0 




(g) 



FIG. 2. (a) Quasienergies e' = e - eqj,(C1) and (g) mean energies 
E of Floquet states as a function on the driving frequency f2 and the 

linear momentum fik. In the plot A p = fifl c /2 = ^A^+// 2 . We focus 
on the solutions that are affected by the ideal rotating SOC driving 
in the small k limit for \i <k A , given in Eq.(|47). For the sake of 
clarity, surface plots show the lower energy bands only (see E3(k) 
and E 3 (k) in the text): upper bands follow by symmetry {s 4 (k) = 
-E 3 (k) and E 4 (k) = -E 3 (k)). In panels (b), (c), (d), (e) and (f), the 
driving frequency Q is fixed and both bands are shown. At the critical 
frequency — panel (d) — the quasienergies form a Dirac-cone while 
the mean energies show a similar behavior except for a discontinuity 
at k = 0. lust away from Cl c , sub- and super-critical regimes can 
not be distinguished from the quasienergy bands alone. In contrast, 
thanks to the discontinuity developed at the critical point, the mean- 
energy bands show remarkable features distinguishing one regime 
from another. Such a contrast is clearly seen already in the surface 
plots of panels (a) and (g). 

of the operator obtained by enforcing Ez = in Hq. Similarly, 
two other indicators of the band inversion can be constructed 
here: the expectation value of Hq evaluated at Ao = 0, 



E a ~ A, 



dEa 



and the expectation value of Hq evaluated at \i — 0, 

dE a 



E a -fi 



djji 



(51) 



(52) 



Such expectation-value indicators, 1^ (here f = {E z ,fi, Ag}), 
can be constructed from the eigenvalues E a with great gener- 
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ality to deal with band inversions. Physically, this is not sur- 
prising since the complements of namely, E a - 1% are the 
expectation values in the eigenstate of operators appearing in 
the Hamiltonian. As discussed in Sec|H|here the k « eigen- 
states of the bands involved in the topological phase transition 
change abruptly when passing through the critical condition, 
this justifies the observed behavior. Here, we do not discuss 
these indicators further as they are beyond the scope of this 
work. 



4. Other forms ofSOC driving beyond IR case 

We have shown that the topological system of Ref. |6]is di- 
rectly related to Floquet system composed of a non-magnetic 
superconducting wire subject to a rotating SOC axis. The 
studied time-dependent Hamiltonian satisfies: (i) Fourier de- 
composition has only zero frequency and +Q components, (ii) 
the zero frequency component of the Hamiltonian has no SOC 
contribution, (iii) the ideal SOC rotation introduces a spin 
symmetry around the 3-axis, and (iv) the choice of one out 
of two possible rotational senses for the SOC axis determines 
the breaking of time-reversal symmetry (needed for emer- 
gence of Majorana fermions at the boundary of the topologi- 
cal region).B3Here, we briefly discuss the Floquet quasiener- 
gies and mean energies in systems satisfying only some of the 
above conditions. Notice, however, that a more rigorous study 
of the situations below would require analyzing each of their 
effective Hamiltonians (see Eq.(|4T)). 

In Fig. [3] we present a sketch of the ideal rotating SOC, 
panel (a), together with two alternative types of SOC drivings, 
panels (b) and (c). Figure |3jb) depicts a linear (L) driving 
Al(?X sucn that 

A L (r) • <r = A IR (f) • er + o-[ (Air(0 ■ <t)o-i, 
=> A L (f) = (25 cos (Qt) , 0, 0) . (53) 



(a) A ir(0 



(b) A L (0 




(c) a rW 



FIG. 3. Sketch of three different cases for the periodic SOC driving 
A(f) of Eq. (13). (a) Ideally rotating (IR) driving A m (f) of Eq. (14) : 
the SOC axis rotates uniformly in the 1-2 plane with constant cou- 
pling strength, (b) Linear (L) driving AlU) of Eq. l |53} : the orien- 
tation of the SOC is fixed while the coupling strength oscillates har- 
monically, (c) Ramp (R) rotation driving A R (r) of Eq. l |55[ >: the SOC 
axis rotates anharmonically in the 1-2 plane with varying coupling 
strength driven by symmetrical triangular waves 



This represents the superposition of clockwise and counter- 
clockwise IR rotations with the same amplitude a. This driv- 
ing does not fulfill property (iii) of the IR case since the SOC 
strength oscillates while the axis direction is preserved. The 
Fourier amplitudes are 



^(±i) _ £(±i) 



e.L 



*h,L 



a a 
= 2-p x ai =2-p x (o- + +o-_). (54) 
n n 



This means that the reduced Floquet operator for the L driv- 
ing mixes not only states |0, +, k {) and |1, — , k f ), as in the IR 
case of Eq.(|44l), but also states |0, +,k f> and \l,-,k I). This 
extra mixing is a consequence of the violation of property (iv), 
leading to the double-degeneracy of the full quasienergy spec- 
trum, as we shall see below. 

Figure |3jc) depicts a different kind of rotation defined by 



A R (?) = (acos R (Qf) , -asin R (Qt) , 0) . 



(55) 



which is based on a 27r-periodic ramp (R) or symmetrical tri- 
angular wave function 



sin R (;c) = | " ( 4 1) " (x - mn) , for - | < x - mn < | 



- Z 



(-1)" 



n- 



sin(nx), 



n>0,odd 

cosr(x) = sinR(jc + 7t/Z). 



(56a) 



(56b) 



With this choice — notice that cos R («7r) = (— 1)"tt 2 /8 — , the 
contributions at frequencies +Q. are identical to those for IR 



case with Fourier amplitudes ^^,1=^ given in Eq. (|27 1 
Higher frequency Fourier components are 



e,R h,R 



(-1)" 



(57) 



with n are positive, odd integers. The R driving, hence, vio- 
lates properties (i) and (iii) of the IR case. 

So far we have considered three different types of driving 
A(f), without any static SOC component A (()) = aro/Mo (with 
strength ao and spin axis along direction ho), ho can be either 
parallel or perpendicular to the plane defined by A(f). In the 
following, we focus on the case of dominating A(f) such that 
ao ~ Q.2& (nothing particular is expected for ao » a since 
the system behaves, essentially, as a non-topological wire with 
static SOC and vanishing magnetic field). In Table[l]we define 
the labels used to identify 9 different situations addressed nu- 



merically with the method discussed in Sec. Ill A| Notice that 



drivings L and R, together with most cases with nonvanishing 
A (0) , cannot be easily mapped to a static model as in the IR 
case. 

The results are presented in Fig. |4] At the critical point 
Q = Q c , the Floquet-band gaps close at s — eo,o(^c) = Sf2 c /2 
for k — > 0. For L driving, as in the IR case, Floquet-band gap 
(at small k) closes for the critical frequency forming Dirac- 
like cones. However, the bands are double degenerated and 
protected edge states are not favored: in a finite system two 
edge states would actually appear at each edge, but subject to 
a strong mixing by local perturbations. Here even in absence 
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FIG. 4. Summary of results for the situations listed in Table[T]with /j = 1.35A , TiCl c « 3.36A , and E so = a 2 m* /(2Jf) = 0.15A . Whenever present 
(lower panels), the static SOC component is set ao = a/5. Labels subindeces "<", "c" and ">" stand for fl = 0.9fJ c , f2 c , and l.lfl c , respectively. 
Floquet quasienergies are shown from to TlQ c . Floquet bands are symmetrical with respect to s = HQ/2 (horizontal grey line). We do not show 
the —k quasienergies and mean energies whenever they are identical to the +k ones. This is not the case for finite static SOC orthogonal to the 
periodic SOC driving. For vanishing static SOC (upper panels), solid black lines depict IR-quasienergy Floquet solutions — in particular, the 
bands associated to the topological transition [quasienergies from Eq.l|33|> with n = generated by the energies of the static problem £+ _(£) and 
E- + (k) given in Eq.(|3 ![>]■ For L driving, the solutions with ao =0 are doubly degenerated. For R driving, higher-frequency driving components 
generate new transitions. Floquet transitions introduce avoided crossings in the quasienergy spectrum and strong discontinuities in the mean 
energiespEHSee text for further discussion. 



of perturbations the localized states do not appear because 
those states mix each other. For completeness, in Figj4]we 
also present results for the L driving coexisting with a static 
SOC where the referred double degeneracy is lifted. However, 
we have checked that no edge states can be found in the asso- 



ciated finite system (Sec IIIB i. This is something expected as 
we find two Floquet bands with finite slope closing the gap at 
the critical frequency. 

For R driving, in the vicinity of k = the Floquet bands are 
similar to those for the IR case, except at the critical frequency 
where a second pair of bands with flat dispersion also close. 
For larger k, the effect of higher harmonics modifies signifi- 
cantly the shape of quasienergy bands. The avoided crossings 
in quasienergy produce abrupt changes in the mean energies, 
the spikes in mean energ ies gro w the larger the amplitude of 
the excitation.EElln Sec. Ill B we show that related finite sys- 
tems can hosts Majorana solutions at the edges, but presenting 
some signatures of the higher harmonics. Considering that the 



amplitudes of the R-driving Fourier components (Eq.(|57|)) de- 
cay very fast, these results suggest that Majorana fermions at 
the edges might be absent for arbitrary SOC drivings with sig- 
nificant contributions from higher harmonics. 

Finally, we briefly discuss the effect of a static SOC com- 
ponent with results shown in Fig. [4] For both IR and R driv- 
ing, a A (0) orthogonal to the plane defined by the driven SOC 
axis produces a breaking of +k symmetry. Despite symmetry 
breaking, the main features discussed for vanishing A (0 ^ are 
maintained. This is related to the fact that the IR-mapping 
U^(t) leads to static model of Ref. |6] except for a small SOC 
component along the magnetic field axis. 

Besides, for parallel IR driving, closing Floquet bands at 
critical frequency show a tendency to reduce their mean en- 
ergy. Moreover, the presence of Majorana states in cor- 
responding finite systems is discouraged (see Sec. IIIB I 



whereas the orthogonal IR driving shows bounded MF undis- 
turbed as those for oq =0. 



B. Finite systems: Floquet Majorana Fermions 



TABLE I. Labels used to identify different SOC drivings A(f) in the 
presence of a static SOC component A <0) . 



A (0) 


A(r) 


Air(0 


A R (f) 


A L « 


A (0) 


= 


IR 


R 


L 


A (0) 


II A(0 


IR" 


Ril 


L» 


A (0) 


J.A(f) 


IR X 


R x 


L x 



The Floquet Majorana fermions expected for this system 
have the form given in Eq. ([36| with quasienergy efmf = 
£o,n(£2) = (n + l/2)ftfl Its linear dependence on the driving 
frequency determines that the mean energy of FMFs is zero, 
see Eq. < [23] >. In the following, we apply a numerical method 
looking for solutions satisfying s « eo,«(fi) and E « in fi- 
nite systems subject to a variety of drivings. Our numerical 
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approach solves the eigenvalue problem for the evolution op- 
erator over a driving period, obtaining the quasienergies based 
on Eq.(fl~8]>. As the evolution operator is evaluated from t — 
to t = T, its eigenvectors are the instantaneous form of the 
Floquet states at t = 0. 

We compute the evolution operator by dividing the period 
T in 120 time slices during which the driving is assumed to 
be constant. Due to the finite sample size, the system must be 
treated in real space: we use the customary finite-differences 
passage from the continuos to a tight-binding model with N 
sites. The wavefunction at each site is represented by a Nambu 
spinor and the pairing potential enters as an on-site coupling 
between the electron and hole sectors. The evolution opera- 
tor becomes a AN x 4-N matrix, Ujit). The lattice spacing oq 
is chosen sufficiently small to turn the kinetic term [second 
order spatial derivative, now a hopping fh = ft 2 / (2m* a?) be- 
tween neighbor lattice sites] into the largest energy scale in the 
system. The SOC [proportional to first-order spatial deriva- 
tives] leads to nontrivial spin-hopping terms with amplitudes 
given by the components of the vector (A(?) + A (0) ) /(2an). As 
usual, the tight-binding model has the flexibility to deal with 
any site-dependent perturbation as local disorder or inhomo- 
geneous driving, both discussed in Figs. [5] and [6] 

We start by studying a clean, finite sample subject to IR 
driving in the supercritical regime Q. > Q c . Bound states at 
the edges are expected as a result of change in topology, this 
is because the vacuum is in the trivial phase: it can be inter- 
preted as having a large and negative chemical potential and 
therefore the critical frequency becomes larger than the driv- 

.5 



a) we show the re- 



ing frequency, HQ. < 2 ^u^-t-A^. In Fig. 
suiting quasienergies and mean energies. We sort the solutions 
according to growing quasienergies (left panel) or to growing 
mean energies (right panel), finding two solutions which are 
close to the Majorana conditions 
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FIG. 5. Floquet Majorana states in a finite sample subject to the 
IR excitation in the Q. > Q c phase with E so = a 2 m*/(2% 2 ) = 0.5Ao. 
Reference quantities are A and the wavelength A A = 2n/k A , with 
k A = V2"i*Ao/fi. The lattice spacing is «o = 0.0225/1a (see text) 
with sample length 300«o = 6-7 5 A A . We choose \i = 2A such that 
Kl c = 2 V5A , with driving frequency fl = 1.15fS c . (a) Floquet solu- 
tions ordered by growing quasienergy (left) and mean energy (right) 
. We plot e' = e-h£l/2 for < e < TiCl. Two edge states, \a) and \b), 
are found with s as TiCl/2 and mean energy E ss 0. Isolated Floquet 
Majorana fermions at each edge, |<* r ), are extracted from the bond- 
ing and antibonding combinations of \a) and \b) (see text). For the 
two FMF states at t = we plot, as a function of the position x: (b) the 
probability density, \cf>l | 2 = \(f> T (x, t = 0)| 2 ; (c) the spin density along 



0\l(O) = ""2. $ up) = ^±(t)( x ' 1 = m ^ W tne s P^ n density along tr 3 , 
st" = 5* (x, t = 0). The spin density along <X|| (0) = cr l (not shown) is 
zero, (e) We include uncorrelated disorder E, e [-0.75Ao, 0.75Ao]. 
FMF are found with probability density similar to that of clean sys- 
tems shown in (b). FMF densities are plotted as a function of time 
along one driving period in the presence of disorder: (f) \(p T (x, t)\ 2 , 
(g) S^.Ax, t), and (h) 5* (x, t). The in-plane spin component rotates 
while the FMF state conserves its location (see text). 



e=s ,„(Q), £ = 0. 



(58) 



We call these solutions \a) and \b), both found to have prob- 
ability weight at the two edges. Similar features are found in 
Ref. when calculating the eigenstates of finite, static sys- 
tems. These eigenstates are bonding and antibonding com- 
binations of Majorana fermions at the edges produced by fi- 
nite size effects. The eigenenergies +5E small but different 
from zero. To get rid of the mixing we should undo the 
(anti)bonding as 53 
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i>l(t=0))=^=(\a)±\b}). 



(59) 



In this way (see Fig|5jb)) we obtain two independent Flo- 
quet solutions, \<p+(t = 0)j and \cp^(t-0)j, one on each sam- 
ple's edge. These are the solutions actually expected on long 
samples where the overlap between localized states is expo- 
nentially small. 

For studying the structure of the these solutions we first in- 
troduce the electron-hole identity and Pauli matrices r ; - with 



j = 0, 1,2, 3. In this way, the operators in the 4-dimensional 
Nambu space can be written as a combination of the operators 



TjCT; = 



T T j,220~i J 



(60) 



where Tj^g (y, 6 — 1,2) are the components of the r, matrix 
such that Tjy g o-i are 2 x 2 blocks operating on spin space. The 
components [1,2,3,4] in the Nambu space are associated to 
to" = [++,H — , — h, — }, respectively, where r and cr refer to 
the ± 1 eigenvalues of T3 and <r 3 . Notice that the combination 
T e = (to + Tj)/2 [t/j = (to - Ty)/2] has one single nonzero ele- 
ment, Te n = 1 [i"/?,22 = 1], and it is useful for writing operators 
associated to the electron [hole] sector only. 

For the time-dependent solution |^(f)) = e~s E '|0 r (/)) we 
define the instantaneous probability density \<p T (x, f)| 2 as the 
sum of the electron and hole probability densities of the as- 
sociated Floquet state, p e i h (x, t). These are obtained when 
computing the expectation value of the identity operator 1 = 
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f dx\x) {x\ {T e +T h )a-Q\ 

dx\pf(x,t)+p((x,t)j, 
pf /h (x, OeJ( </> T (f)\ xro) <T|r e// ,|r') (jctV \<p T (t) ) . (61) 

Similarly, related spin densities also evolve in time. The to- 
tal spin density, S? t( Jx,f), is the sum of the electron and 
holes probability densities in the associated Floquet state: 
Sf e i h (x, t). These are obtained by computing the expectation 

value of S i = J dx\x) (x\ (T e +Th)o~i, with i= 1,2,3 (H/2 factors 
skipped for simplicity): 

dx[s( e {xj) + sf h {xj)y (62) 

s tih^ f ) =Yj ( ^ T{t) \ xto ) to T i^/A cr 'i tr '' r '> ( xrV k r w ) • 

CT.tx' ,t,t' 

Notice that the densities defined in Eq.(|6T]l and Eq.(|62]> are 
properties of the physical state because the quasienergy phase 
factor cancels out. Moreover, due to such cancelation the re- 
sulting densities are the same irrespective of which shifted 
version of the Floquet state (see Eq.(|T9)>) is taken. 

As mentioned above, the states |</>±(f)) studied in this sec- 
tion [constructed with Eq.([59]>] satisfy (numerically) the Ma- 
jorana condition of Eq.(|58]l. From hereon, we use |0 r (?)) to 
refer to any of these two states. In all presented situations 
(Figs, [i] and ^ we find \(f> T (x, t)\ 2 /2 = pf(x,0) = pf (jc.O) 
and Sf tot (x,0) = for all i = 1, 2 and 3, namely, sf(x,0) = 

—S?,(x, 0). Notice this differs from the regular type of quasi- 
particle states arising in an s-wave spin-degenerated BdG 
Hamiltonian; the latter cannot have vanishing ST. (x, 0) for 
all three spin directions simultaneously. However, the result 
coincides with what is expected for a Majorana fermion as its 
own antiparticle (see Nambu-space representation of Eq.(|36])). 
In what follows, we focus in the spin densities for electrons 

Sf (x,t) = 5 1 . (x,i) as it is meaningless to discuss the total 
spin of a Majorana solution. 

In Figj5] for both clean and disordered systems subject to 
IR driving in the Q. > Q c phase, we present the total prob- 
ability density \(f> T (x, t)\ 2 and the electron spin densities. In 
Fig j5jd) we see that the states at different edges have the same 
S 3 (x, 0) profile. Moreover, these profiles do not change with 
time as shown in Fig|5Jf). This is equivalent to what found for 
the electron spin density along the magnetic-field axis for MF 
in the static system of Ref. 6 (where the larger the SOC in rela- 
tion to the Zeeman energy, the smaller this component^). In 
our non-magnetic system, the associated direction is normal 
to the plane defined by the IR driving, see Eq.(fT4|i. 

Additionally, in Fig|5jc) we see that the FMF solutions have 
opposite electron spin in-plane components which change as 
a function of time as shown in Fig. |5|g). We have defined 
the parallel, (r\\( t ), and the perpendicular, cr ± (f), spin directions 



with respect to the instantaneous SOC axis Ai R (f): 

(T\\( t ) = o~\ cos(Qf) - cr 2 sin(Qf), (63a) 
o\l(0 = o"2 cos(Qf) + o-\ sin(Of)- (63b) 

We find that the spin density StJx,t) is zero (not shown) 

while S^(f)(x, t) in Figjijd) remains constant. This indicates 
that the in-plane electron spin components rotate as a function 
of time. This rotation is already expected from Eq. |36]l as a 
mapping of the spin properties of MFs present in the static sys- 
tem of Ref. 6 , which show a finite electron spin density nor- 
mal to both the SOC axis and the magnetic-field axis. Such 
spin component assumes opposite signs for MFs located at 
opposite edges of the sampled in agreement with our results 
entirely obtained within the periodically driven system. 

Figure |5je) presents one realization of uncorrected disor- 
der incorporated in all remaining results of this section, in- 
cluding time-evolved properties of Floquet states depicted in 
Figs.|5|f),|5jg) and|5Jh). The disorder is implemented as tight- 
binding onsite energy fluctuations E, e [-0.75Ao, 0.75Aq] 
chosen from an uniform distribution. Its characteristic energy 
is comprable to A () and larger than the SOC-driving energy 
scale excitation, £ so . We have deliberately used short-range 
disorder in order to avoid that long sections of homogeneous 
disturbances drive parts of the system outside of the topologi- 
cal phase to be characterized. The robustness of the results to 
disorder is a signature of their topological origin. 

Figure [6] shows several situations alternative to the IR driv- 
ing including disorder. There we plot the probability densities 
for the states |<^ r (f)) at f = for a supercritical driving fre- 
quency O = 1.15Q C . The Fig(6ja) presents results for R driv- 
ing. Localized states at edges have different profiles due to the 
presence of higher harmonics, evolving as a function of time 
in contrast to IR driving (not shown). 

Systems subject to L driving (not shown) do not present 
Floquet states satisfying the Majorana conditions of Eq. ( f58] l. 
This is due to the double degeneracy of Floquet states, in- 
consistent with topologically protected unpaired edge states. 
Similarly, the presence of partial L driving tends to delocal- 
ize Majorana states. In Fig j6jb) we show a system subject to 
hybrid IR/L drivings on the left and right of the sample, re- 
spectively. We find that one of the states is localized on left 
edge while the partner is delocalized all over the sample. The 
presence of IR driving introduces Majorana-like solutions but 
only one of them is localized at the edge because of the L 
driving. 

Figures|6|c) and|6|d) show results for IR driving in the pres- 
ence of an small static SOC component, both parallel and or- 
thogonal to the plane defined by the driving, respectively. In 
either case, localized Majorana solutions similar to those for 
the pure IR case are found at the sample edges. The orthogo- 
nal configuration is less affected by the static SOC, presenting 
a highly symmetric state distribution in contrast to the paral- 
lel configuration. This is consistent with the behavior of the 
mean energies (cases IR" and IR -1 in FigQ for extended solu- 
tions with small k. 
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FIG. 6. Floquet Majorana states in finite samples subject to different types of SOC drivings and/or static SOC. All samples include disorder 
as the one given in Fig.[5je). Unless otherwise stated the parameters cig, fi, Q. c , f2 = 1.15fl c and E so are those given in Fig. [5] As in Fig(5] 
decoupled solutions are obtained by combining Floquet states |a> and \b) with EK.hQ.j2 and £*0. In all cases we plot the probability density 
l^jl 2 = |0 r (x, t = 0)| 2 of the two (or more) FMF candidate states. We rescale some curves (x5) to improve visibility, (a) Full sample with 
R driving, (b) Left (/) section of the sample is subject to IR driving whereas the right section (r) is L-driven. (c)-(d) IR-driven sample in 
the presence of a static SOC with ao = 0.2a either parallel (c) or perpendicular (d) to the plane of the IR driving, (e) IR-driven sample with 
inhomogeneous chemical potential on the left and right sections, // = 2A = n' 12. The driving frequency is 1 . 150J. but smaller than the £l r c . (f) 
IR-driven sample on the left section in contact with no SOC on the right, (g) IR-driven sample subject to different driving frequencies on left 
and right sections, Q; = 2£l r = 1.15£2 C : the left (right) section is in the supercritical (subcritical) regime, (h) Sample subject to counterrotating 
IR drivings with Q.= 1.15Q C , clockwise (counterclokwise) on the left (right) section of the sample. Four Floquet states arise in this case. 



In Fig[6je) we present the results for IR driving with in- 
homogeneous chemical potential. The driving frequency is 
chosen to supercritical on the left and subcritical on the right. 
Majorana states are found at the edges of the left section. This 
is in full agreement with what expected by virtue of the map- 
ping to the static topological system of Ref. |6] 

In Fig j6ff) we show the Floquet Majorana states for the case 
in which the left section is subject to IR driving (Q > D. c ) 
while the right part is free of any SOC. We find two Majo- 
rana states at the edges of the left section, this time with one 
delocalized solution penetrating into the right section. This is 
expected, indeed, since the section with vanishing SOC maps, 
using the time-dependent spin rotation for the whole system, 
into a static superconductor in which the Zeeman energy is 
larger than Ao, leading to E = solutions that are not gapped 
because the SOC is zero. 

In Fig. |6|g) we have investigated a situation in which the 
sample is subject to inhomogeneous IR-driving frequencies: 
f2 r on the right section and Q 1 = 2£2 r on the left one. We 
choose a supercritical fi 1 = 1.15Q C such that Q 1 is subcriti- 
cal. We avoid dealing with incommensurable frequencies so 
that our choice allows for the application of the Floquet the- 
orem. The results are similar to the case of Fig. [6]T), with a 
localized Majorana solution on the left and a delocalized state 
penetrating into the right section. 

Finally, we obtained the Floquet solutions for counterro- 
tating IR drivings: clockwise rotation on the left section and 
counterclockwise on the right one, see Fig. |6|h). The situa- 
tion results interesting since both sides satisfy the supercritical 
condition Q. > £2 C . As in all other cases, we have chosen the 



position of the left/right interface away from the geometrical 
center of the sample to avoid solutions influenced by symme- 
try. We find two pairs of states satisfying the Majorana con- 
dition. The finite-size offsets to Eq.(|58]>, ±6s and +6E, are 
different for the two pairs. For each pair we find a localized 
FMF at one edge and a delocalized solution over the opposite 
side of the sample. 

The results discussed along this section, obtained in the 
presence of disorder, confirms the existence of bounded Flo- 
quet states at the edges systems subject to IR driving. They 
are Floquet Majorana fermions arising in absence of magnetic 
fields when the driving frequency is sufficiently large. When 
the samples are enlarged, the local probability density of the 
delocalized FMF solutions decrease [see Figs.|6jb),[6jf),[6jg), 
and[6|h)] whereas the localized solutions remain unaffected. 

IV. CONCLUSIONS 

Starting from a well-known quantum system hosting Ma- 
jorana fermions for sufficiently large magnetic fields, we de- 
rived a non-magnetic scheme sharing such properties thanks 
to periodic driving. Our proposal disregards Zeeman cou- 
pling: its role is played by the driving frequency of a rotating 
spin-orbit coupling axis. An effective breaking of time re- 
versal symmetry is induced by the definite rotation sense (ei- 
ther clockwise or counterclockwise). Otherwise, double de- 
generacy (due to a time-reversal symmetry) would discourage 
the appearance of a single Majorana solution at the sample's 
edges. A experimental realization of this theoretical model 
would require the ability to control the SOC axis as a func- 
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tion of time. This can be envisioned in systems with lateral 
gates in orthogonal directions that modulates the symmetry of 
the quantum well containing a superconducting electron sys- 
tem. Alternatively, in cold-atom quantum wires the Rashba 
coupling could be engineered to change its axis with timePSH 
The proposal incorporates the virtue of being a non-magnetic 
platform for Majo rana ferm ions into the active field of Floquet 
topological matter! 9 ! 18 ! 21 ! 24 " 27 ! 

Our mapping allows analytical treatment in the case of ideal 
rotating (IR) driving. For supercritical driving frequency Q. > 
Q. c , i.e., in the Floquet topological superconducting phase, 
Floquet Majorana fermions are found at the edges of finite 
samples. While in static systems Majorana fermions appear 
at zero energy, the FMFs here have quasienergy hQ.(n + 1/2). 
We point out that these FMFs, as those reported in Ref.|9]hav- 
ing quasienergy linear in O, have zero mean energy. This fact 
can be important on the light of Ref. |38] where Arimondo et 
al. provide experimental and theoretical indications that the 
mean-energy variable — once the memory of the initial condi- 
tion is lost — might play a key role in dictating the occupancy 
of Floquet states. 39 As the chemical potential lies at zero en- 
ergy, this might determine that, at zero temperature, the re- 
ported FMFs are the highest mean-energy occupied Floquet 
states. 

We have also investigated different types of drivings based 
on time-dependent SOC. Rotating SOC drivings are needed 
because of the requirement of time-reversal-symmetry break- 
ing. For non-ideal SOC rotations the numerical results show 
that higher harmonics can degrade the formation of FMFs. 
Similar degradation is found in the case of an IR driving coex- 
isting with a static component of SOC contained in the plane 
of the rotation. We demonstrate the great deal of possibilities 
for FMFs at edges and interfaces by exploring only a handful 
of examples. 

For finite samples, our numerical simulations show that the 
FMFs appear even in the presence of disorder. This is re- 
markable, considering that FMFs in a finite sample are not 
protected by a gap in quasienergy. In the case of IR driving, 
numerical results can be understood analytically by a mapping 
to the static model, noticing that static disorder maps trivially: 
Majorana fermions in those systems are robust to disorder 
smaller than the protecting energy gap. Interestingly, an un- 
derlying "gap", which is not visualized straightforwardly from 
the Floquet quasienergy spectrum, is protecting the FMFs in 
these driven systems. 

Finally, we expect that tunneling probing (with a normal 
lead) of the edge of a sample in the Floquet TSP phase will 
lead to observe two peaks in the differential conductance 
dl/dV at bias voltages V± = +MI/2. Those peaks have dif- 



ferent heights dictated by the ratio between the effective Zee- 
man energy and the spin-orbit coupling strength. The peaks 
are not present in the trivial phase (Q < Q c ) and a background 
with structure is visible in both phases (including pumping 
effects,™ as current for zero bias due to the rotating-SOC). 
These qualitative conclusions are obtained within a picture 
that makes use of the rotating frame transformation [the driven 
part of the sample is mapped to the static topological phase, 

having an unpaired MF at its edge at E — 0] and shifts the 
Fermi energies in the normal probe differently for electrons 
with spin cr z = +: an energy cr z Ml/2. 59 The mentioned shifts 
in the Fermi occupations are introduced to conserve the num- 
ber of electrons in the static lead before and after the transfor- 
mation (notice that the effective magnetic field splits the up 
and down bands) as the ultimate goal is computing transport 
from the lead to other parts of the system.^ 

Certainly, among future directions of study appears a quan- 
titative determination of the mentioned observable signatures 
(and others) to detect the Floquet TSR Results beyo nd the 
IR case would requiere a Floquet-Keldysh approach ] 21 ! 61 ! 62 ! 
which here becomes non-standard because the superconduct- 
ing ID lead has a time-dependent spin-orbit coupling. Other 
line includes developing schemes to use the localized Floquet 
Majorana fermions (the ones reported here or those reported 
in Ref.9) for quantum information purposes: quantum mem- 
ory, braiding, etc. Crucial in all cases appears having a better 
understanding of which of the Floquet states are occupied, a 
problem that falls within the complex subject of statistical me- 
chanics for driven systems! 6 - 3 J 6 ^lHere, we have pointed out that 
interesting intrinsic physical properties are expected if one de- 
velops the ability to manipulate the SOC axis as a function of 
time. Identifying candidate platforms in which such ability 
is achievable would be desirable and it certainly would mo- 
tivate more research on this and other systems, this includes 
situations in absence of superconductivity. 
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